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SEMI-CLASSICAL ANALYSIS OF A RANDOM 
WALK ON A MANIFOLD 

By Gilles Lebeau and Laurent Michel 

Universite de Nice Sophia- Antipolis 

We prove a sharp rate of convergence to stationarity for a natural 
random walk on a compact Riemannian manifold (M,g). The proof 
includes a detailed study of the spectral theory of the associated 
operator. 



1. Introduction. This paper has two main aims. First, we study the 
spectral theory of a Markov chain associated to a natural "ball walk" on 
a compact, connected Riemannian manifold. From x, the walk moves to a 
uniformly chosen point in a ball of radius h around x. Here h is a small 
parameter. We prove a precise Weyl-type estimate on the number of eigen- 
values close to 1, and convergence of the spectrum near 1 (when h — > 0) to 
the Laplace-Beltrami spectrum. This walk does not have, in general, the 
Riemannian area distribution as stationary distribution. The second aim is 
to analyse the Metropolis algorithm as a way to achieve uniformity. Sharp 
rates of convergence for the Metropolized chain are given. In the Appendix, 
we prove that under appropriate scaling, the modified Metropolis chain con- 
verges to the Brownian motion. 

Let (M,g) be a smooth, compact, connected Riemannian manifold of 
dimension d, equipped with its canonical volume form d g x. Let d g (x,y) be 
the Riemannian distance on M x M. For x £ M and h > 0, let B(x,h) = 
{y,d g (x,y) < h} be the ball of radius h centered at x, and let \B(x,h)\ = 
Ib(x h) dgV be its Riemannian volume. For any given h > 0, let T\ x be the 
operator acting on continuous functions on M, 
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We denote by Kh the kernel of Th, which is given by 

(1.2) Kh ^ y)d ^ V= 1{ \^xM daV - 

Obviously, for any x £ M, Kh(x,y) d g y is a probability measure on M, and 
therefore Kh is a Markov kernel. It is associated with the following natural 
random walk on M: if the walk is at x , then it moves to a point y G B(x,h) 
with a probability given by Kh(x,y) d g y. 

Let t Tf l be the transpose operator acting on Borel measures on M, defined 
as usual by ( T^(/i),/) = (fj,,Th(f)}- Let Cd be the volume of the unit ball 
of the Euclidean space Mr. For h small, h~ d \B(x, h)\ is a smooth function 
on M which converges uniformly on M to q when h — > 0. Let be the 
probability measure on M, 

, \ , \B(x,h)\ , 

(1.3) h Z^J?~ " x ' 

where the normalizing constant Zh is such that dvh(M) = 1. Then for /i 
small, di^ is close to d g x/ Vol(M) and Z/j is close to Vol(M). One verifies 
easily that is self-adjoint on the space L 2 (M,dvh), and that t Ti l (dui l ) = 
du h . 

The first goal of this paper is to analyze the spectral theory of the self- 
adjoint operator T/j acting on L 2 (M,dvh). Let us recall some basic facts. 
One has T/j(l) = 1, and by the Markov property, the norm of acting on 
the space L°° is equal to 1; by self-adjointness, the norm of Th acting on the 
space L?~{M,dvh) is equal to 1 and thus the norm of Th acting on the space 
L 2 (M, dvh) is also equal to 1. Observe that for any given h > 0, the operator 
Th is compact. Thus the spectrum of Th, Spec(Th), is a closed subset of 
[—1,1] which is discrete in [—1,1] \ {0} with as accumulation point, and 
each fj, G Spec(Th) \ {0} is an eigenvalue of finite multiplicity. 

We denote by A g the (negative) Laplace-Beltrami operator on (M,g), 
and by = Ao < Ai < A2 < ■ • ■ < X n < ■ ■ ■ the spectrum of the self-adjoint 
operator — A g on L 2 (M, d g x). We will denote by Gd{C) the following function 
of f € R d : 

(1.4) G d (0 = ~ f e^dy. 

Cd J\y\<l 

Up to the factor — , the function Gd is the Fourier transform of the charac- 
teristic function of the unit ball in M. d , and depends only on |£| 2 . We shall 
also use the function r^(s) on [0,oo[ defined by 

(i-5) G d (0 = r d m 2 ). 
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The function Ta is real analytic, |r^(s)| < 1, and lim^oo Td(s) = 0, since 
Gd(0 is the Fourier transform of a compactly supported, real and even L 1 
function of total mass 1. One has near s = 0, 



(1.6) 



i 



+ 0{s 2 ). 



2{d + 2) 

Moreover, there exists 70 < 1 such that r^(s) £ [—70, 1] for all s, and one has 
r^(s) = 1 iff s = 0. To see this point, just observe that if = 1, then 

one has G^(^) = e l9 for some real 9, hence i| y |<i( eJ2/ ^~* e — l)dy = which 
implies y£ — 9 £ 2ix7L for all \y\ < 1, and therefore £ = and 9 £ 27rZ. 

Theorem 1. Let ho > 6e small. There exist 7 < 1 such that for any 
h£]0,ho], one has Spec(Th) C [—7,1], and 1 is a simple eigenvalue of Tf t . 
Let 



1 



(1.7) < • • • < nk+i(h) < fi k (h) <---<vn(h)< no{h) 

be the decreasing sequence of positive eigenvalues ofT^. For any given L > 0, 
there exists C such that for all h £ ]0, ho] and all k < L, one has 

l-fi k (h) X k 



(I- 



<Ch 2 . 



h 2 2{d + 2) 

Let N(a, h) be the number of eigenvalues of Th in the interval [a, 1] . For 
any given 5 £]0,1[, there exist Cg^ independent of h £ ]0, ho], such that the 
following holds true: 

For any r £ [0, (1 — 5)h~ 2 ], N(l — rh 2 ,h) satisfies the Weyl law, 



(1.9) 



N(l-Th 2 ,h)-(2irhy 



r d (K|2)e[i-rhM] 



dx d£ 



<C S a{1 + t 



,(d-l)/2 



where dxd£ is the canonical volume form on the symplectic manifold T*M, 
and \£\ x is the Riemannian length of the co-vector £ at x. In particular, one 
has 

(1.10) N(l-rh 2 ,h) <C 5 ,2{l + r) d/2 . 

Moreover, for any eigenfunction e\ of associated with the eigenvalue 



fj-k(h) £ [5, 1], the following inequality holds true with T k {h) = h 2 (1 
(1-H) \\e h kh°°< 05,3(1 + T k (h)f A \\e h k \\ L2 . 



Let \Ah\ be the positive, bounded, self-adjoint operator on L 2 (M,duf l ) 
defined by 

h 2 



(1.12) 



l-T h 



2(d + 2) 



I A, 
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By (1.8), the two operators |A/J and — A g have almost the same eigenvalues 
in any interval [0, L) independent of h, for h small enough. Our next result 
gives more precise information on the difference of their resolvents for h 
small. Observe that as vector spaces, the two Hilbert spaces L 2 (M,dvh) and 
L 2 (M,d g x) are equal, and that their norms are uniformly in h equivalent. 
We set L 2 = L 2 {M,du h ) = L 2 {M,d g x), \\f\\ L 2 = \\f\\ L ^ M ,d g x/Vol(M)), and if 
A is a bounded operator on L 2 , we denote by ||A||^2 its norm. 

Let F\ and F2 be the two closed subsets of C, F\ = {z, dist(z, spec(— A g )) < 
e}, F2 = {z,Re(z) > A, \ Im(z)\ < eRe(z)} with e > small and A > large. 
Let F = Fx U F 2 and U = C \ F. 

Theorem 2. There exists C,ho> such that for all h £]0,ho], and all 
zeU, 

(1.13) ||(2 - lA^ir 1 -(z + A g ) _1 || L 2 < Ch 2 . 

Remark 1. The error term 0(h 2 ) in the estimate (1.13) is of the same 
type than the error one gets for the difference between discrete and continu- 
ous Laplacian on M. d . However, in our geometric setting, the Ricci curvature 
of M contributes also to the error term (see Lemma 3 below), and to get a 
true discrete Laplacian on the manifold M, one will have to discretize the 
integration process in formula (1.1). Although this is clearly a question of 
practical interest [as well as modification of |A^| to improve the convergence 
in (1.13)], we will not discuss this point in the present paper. 

Observe that when M = (]R/2-7rZ) rf is the flat d-dimensional torus with g 
equal to the Euclidean metric, one has the equality, 

(1.14) T h = T d (-h 2 A g ). 

Thus, in that case, the operators Th and A s have exactly the same eigenvec- 
tors e lkx , and the results of Theorems 1 and 2 can be proved by a simple com- 
putational verification. For a general compact Riemannian manifold (M,g), 
the two operators and A 9 do not commute, and the formula (1.14) is 
untrue. In Section 2, we will use a suitable /i-pseudo-differential calculus in 
order to show that formula (1.14) remains almost true (in a proper sense), 
modulo lower order terms involving the curvature of M. Then, using the 
results of Section 2, we will prove Theorems 1 and 2 in Section 3. Observe 
that the L°° bound (1.11) on the eigenfunctions of | A/J is the exact analogue 
of what one gets from Sobolev inequalities for the eigenfunctions of A g ; in 
particular, this is certainly not optimal, and it will be of interest to know if 
the Sogge estimates (see [14]) for the eigenfunctions of A g are true for the 
eigenfunctions of \Ah\. However, (1.11) will be sufficient for us in the proof 
of Theorem 3. 
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Let us now discuss the second goal of this paper. For any n > 1, let 
K£(x,y) d g y be the kernel of (T/j) n . Then J A K^(x,y) d g y is the proba- 
bility that the random walk associated to starting at x is in the set 
A after n steps of the walk. When n — > oo, the sequence of probabilities 
Kfo(x,y) d g y will converge to the stationary probability dvh{y), but this is 
not quite satisfactory, since on a general manifold M, dvh(y) depends on h. 
Thus, in order to get a Markov chain with the fixed stationary probability 
d[iM = dgxj Vol(M), we modified the kernel Kh{x,y) d g y, according to the 
strategy of the Metropolis algorithm, in the following way. Let 

(1.15) M h (x,dy) = m h (x)5 y = x + K, h (x,y) d g y, 

where the functions and tC^ are defined by 



(1.16) 



m h (x) = l- )C h (x,y)d g y. 
Jm 



Then, Mh(x,dy) is still a Markov kernel, but now, the operator 



(1-17) M h (f)(x)= / f(y)M h (x,dy) 

Jm 

is self-adjoint on the space L 2 (M, d g x), and therefore one has t Mf l {d g x) = d g x 
for all h. Let M£(x,dy) be the kernel of (M^) n . Our purpose is to get an 
estimate uniform with respect to the small parameter h, on the speed of 
convergence, when n — > oo, of the probability MP(x, dy) toward the invariant 
measure d/j,M = d g xj Vol(M). Let us recall that if p, q are two probabilities, 
their total variation distance is defined by 

lb — Q'IItv = sup \p(A)-q(A)\, 
A 

where the sup is over all Borel sets A. The following theorem tells us that 
this speed of convergence is estimated for h small, as expected, by the first 
nonzero eigenvalue Ai of the Laplace-Beltrami operator- A g . 

Theorem 3. Let ho > small. There exists A such that for all h G ]0, ho] 
the following holds true: 

e -y>(h)nh* < 2 gup \\MZ(x,dy) - d m \Wv, 

(1.18) 

sup \\M%(x,dy) - ^a/||tv < Ae~^ hnh2 for all n. 

Here 7(/i),7'(/i) are two positive functions such that j(h) ~ "j'{h) ~ 2 (d+2) 
when h — > 0. 
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Of course, the analogue of this result is also valid if one replaces Mh 
by T/i and dfiM by dvh, with a simple proof. Theorem 3 will be proved 
in Section 4. We will verify that Mh is a sufficiently small perturbation of 
7)i, and, in particular, that estimates (1.11) and (1.10) remains true for 
its eigenfunctions. Finally, in Theorem 4 of the Appendix, we will answer 
a question of one of the referees of the paper, about the convergence of 
the Metropolis chain to the Brownian motion on the Riemannian manifold 
(M,g). 

Perhaps the main contribution of this paper is the introduction of micro- 
local analysis as a tool for analyzing rates of convergence for Markov chains. 
These result in a fairly general picture; the top of the spectrum of the 
Metropolis chain converges to a Laplace spectrum. Because of the hold- 
ing, the Metropolis chain has a continuous spectrum but this is bound from 
±1 and does not enter the final result. This picture was found in a simple 
case in [4] and for the Metropolis algorithm in Lipschitz domains, including 
the random placement of N hard discs in the unit square, in [5]. The present 
paper shows that the picture holds fairly generally. Throughout this paper, 
we will use basic techniques in semi-classical analysis, for which we refer to 
[13] and [7]. 

For an introduction to the well-developed area of probability theory on 
Riemannian manifolds we refer to [11]. For the analysis of the Metropolis 
algorithm, we refer to [6] and references therein. There are also emerging 
applications to statistics on Riemannian manifolds (see [1-3, 10] for examples 
and references). All of these applications lead to the problem of drawing 
random samples from the uniform distribution. This topic has not been 
widely addressed. Some algorithms are suggested in [3]. The present paper 
is a contribution to a rigorous treatment, giving reasonably sharp bounds 
on rates of convergence. 

2. The symbolic calculus of T^. We first recall some basic facts on the 
classical /i-pseudo-differential calculus. For m € K, let S m the set of functions 
a(x,£, h) smooth in (x,£) £ M 2rf , with parameter h € ]0, 1] such that for any 
a, /3, there exists C a ^ such that for all £ M 2d and all h £]0, 1] one has 

(2.1) \d%dj!a(x,{;,h)\<c a ,p(i + \t;\r-M. 

For a £ S m , we denote by Op{a) the /i-pseudo-differential operator acting 
on the Schwartz space 5(M d ), 

(2.2) Op(a)(f)(x) = (2nh)- d J e^-y^ h a(x,^h)f(y)dyd^. 

Let us recall that for a £ S°, the operator Op(a) is uniformly bounded in h on 
the space L 2 (R d ), and that for a £ S m , b G S k , one has Op(a)Op(b) = Op(c) 
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where c = a$b G S m+k is given by the oscillatory integral 

(2.3) c(x, C /i) = (2Tvh)~ d [ e- ize/h a(x, £ + 6, h)b(x + z, £, /i) cZ0, 



and admits the asymptotic expansion 



c(x,£,/i) = 2^ -^—^d%a(x,£,h)d%b(x,£,h) 



\a\<N 



(2.4) 

+ ^r w (x,e,/ i ), r N eS m + l ~ N . 

The subset of S m is the set of a(x,£,h) G S"™ such that there exists a 
sequence a n (x,£) G S m ~ n ,n > 0, such that for all AT, one has 



(2.5) a(x,€,h)= ]T (h/i) n a n (x,£) + h N r N (x,ti,h), r n GS m ~ 

0<n<N 



From (2.4), one has a$b G S™ +k for a G 5™ and 6 G S k cl . 

Let (M,g) be a compact smooth Riemannian manifold, and let ej(x) G 
C°°(M),j > 0, be an orthonormal basis in L 2 (M, d s x) of real eigenvectors of 
— A ff with —A g ej = XjCj. For any distribution / G T>'(M), the Fourier coeffi- 
cients of / are defined by fj = J fej d g x and one has f(x) = ^ . fjej (x) where 

the series is convergent in V'(M). For s G R, let £F(M) = (1 - A S )~ S / 2 L 2 (M, 
d g cc) be the usual Sobolev space on M. For / G V(M) one has / G H S (M) 
iff II /II iff* (M) = Ylj(^ + ^j) s \fj\ 2 < 00 • We shall also use the semi-classical H s 
norms defined by 

(2.6) ii/iiL=E( 1+ ^) s i/;I 2 - 

3 

A family of operators Rh, h G]0,1], acting on the space of distributions 
V'{M) is said to be smoothing iff for any s,t, N, Rh maps H S (M) in i7*(M) 
and there exists C s tN such that for all h G ]0, 1] one has 

(2.7) \\Rh(f)\\w(M)<C s , t ,Nh N \\R h (f)\\ ff"(M)- 

A family of operators Ah, /iG]0,l] acting on the space of distributions 
T)'(M), belongs to the set £™ of classical /i-pseudo-differential operators 
of order m, iff for any xq G M, there exists an open chart U centered at 
xq and two functions C^(U) equal to 1 near xq with ip equal to 1 

near the support of <p such that A^ip = ipA^f + Rh, with Rh smoothing and 
there exists a ~ S n >o(^A) na n( x ''0 £ S™, such that in the local chart U, 
one has ijjAh^p = Op (a). The principal symbol of Ah, oo{Ah){x,£), is by def- 
inition the first term cto(x,£) in the asymptotic expansion of a(x,£,h). It is 
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a well-defined function on T*M, and for any smooth function ip G C°°(M), 
one has 

(2.8) e ^( x )/h Ah{e Mx)/ h) = ao{Ah){XjMx)) + (D (h). 

Then £ c i = \J m £™ is the algebra of classical /i-pseudo-differential opera- 
tors on M. For A h G £% and B h G one has A h B h G £"/ +fc , ^(Aft^) = 
a (A h )a (B h ) and the commutator B h ] = A h B h - B h A h satisfies [A h , B h ] G 
hgrn+k-l ; cr ( i B^]) = {(T (A ft ), ct (%)} where {/, g } is the P oisson bracket. 
Moreover, for any A^ G £™, one has A* h G £™, <7o(^4/j) = co(^U)> and for any 
s£K, there exist C s independent of h G ]0, 1] such that 

(2.9) \\A h f\\ h>s . m <C s \\f\\ h)S Vf£H s (M). 

Let us recall that for any $ G Cg°([0, oo[), the operator <£(— h 2 A g ) defined 

by 

(2.10) <S>(-h 2 A g )(f) = Y,*(h 2 \ 3 )fMx) 

j 

belongs to £~i°° = C\ m ^d-> an< ^ ^ s principal symbol is equal to 
(2.H) a (H-h 2 A g )) = m\l), 

where is the Riemannian length of the co- vector £ at x. For a proof of 
this fact, we refer to [7]. 

Definition 1. A family of operators C^, h G]0, 1], acting on the space 
of distributions V'(M), belongs to the class £^ if and only if Ch is bounded 
uniformly in h on L 2 (M) and for any $o £ Co°([()> co^i one has 

(2.12) $ (-h 2 A g )C h and C h & (-h 2 A g ) belongs to £~ l °°. 

Let Td,h be the operator Y^h = Fd(— h 2 A g ), so that 

(2-13) r^C/Xx) = ^r d (/ i 2 A J )/ i e i (x). 

i 

Since $oFd G Cq°([0,oo[), one has obviously T^^h G £^. 

Let U C M be an open chart with local coordinates x = (xi, . . . , x^) G M d . 
Then for x G U and r > small, the geodesic ball of radius r centered at x 
is given by 

(2.14) B(x, r) = jx + u, kij(x, u)uiUj < r 2 |, 
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where (kij(x,u)) is a smooth and symmetric matrix in (x,u) such that 
kij(x,0) = gij(x). For any function / compactly supported in U and h 
small, Tfrf is supported in U and given in these local coordinates by 

(2.15) T h f(x)= / / f(x + u)Jdet(g(x + u))du. 

\o[x, n)\ Jt uk (x,u)u<h 2 v 

Using the new integration variable hv = w = k l l 2 {x,u)u in (2.15), we get 

h d f 

(2.16) T h f(x) = — — — / f(x + hm(x,hv)v)p(x,hv)dv, 

\B{x,h)\ 

where m(x, w) is the smooth, symmetric and positive matrix, such that near 
u = one has w = k l l 2 (x, u)u <^ u = m(x, w)w, so m(x, 0) = g~ l l 2 (x), and 
p(x, w) = y/det(g(x + u)) \ det is smooth in (x, w) and p(x, 0) = 1. 

Lemma 1. For ho > small and any k, is a bounded operator on 
C k (M) uniformly in /ig]0, ho]. Moreover, there exists C independent of h 
such that, with |A^| defined in (1.12), one has for all f G C 2 (M), 

(2.17) \\\A h \f\\ L ~<C\\f\\ C 2. 

Proof. The first assertion is obvious from (2.16) since ^ is a smooth 
function of x,h € [0, /to]- From (2.16) and the Taylor formula f(x + y) = 
f(x) + Vf(x)y + C(y 2 ||/|| C 2), one gets easily that (2.17) holds true. □ 

In the above open chart U, we define the symbol of T/j, cr(Th) by 

(2.18) a(T h )(x,C,h) = e- ix ^ h T h (e ix ^ h ). 

For any compact set K C £/, there exists /i^ > such that a(Th)(x,^,h) is 
well defined for x £ K, £ € M d and /i e ]0, ^] . From (2.18), one has 

(2.19) <j(T h )(x,£,h)= J 1 * [ e^- m ^ hv >p(x,hv)dv, 

\B{x,h)\ J|„|<! 

and therefore, for any a,/3, there exists C Q>y g independent of /i such that 

(2.20) |^af < 7(T h )( a ; ) e,/i)|<C a|j9 (l + |£|) |a| . 
Observe also that, since m(x, 0) = g ,_1 / 2 (x) and p(x,0) = 1, one has 
(2-21) a(T h )(x,H,0) = T d ml). 

Lemma 2. Zei /io small. For h E]0,ho], the operator belongs to the 
class £® t . 
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Proof. Let M = [j k Uk be a finite covering of M by local charts Uk, 
and 1 = Ylk^Pk{x) a partition of unity with tp^ G Cg° (£/&). Let £ Cg°(?7fc) 
equal to 1 near the support of <pf~- Then for h small enough, one has 

(2.22) T h (f)(x)=J2l>kTh(<Pkf)(x). 

k 

Let T/j fc = ip^Th^Pk'i we reduce to show that for any k, T^k G £®,. Let 3>o £ 
Cq°[0, oo[; there exists tp G C^^Uk) and a compact set K C such that 
Vfc^o( — h 2 A g ) = Op{a)ip + Rh with a(x,£,h) G S 1 ^ 00 with support in £ G if, 
and Rh smoothing. By Lemma 1, T^Rh is smoothing, and thus we are reduce 
to show that in the local chart U^, one has Tj l Op{a) G ■ From (2.2) and 
(2.16), one has 

T h Op(a)(f)(x) = (2Kh)~ d j e^-y^ h b(x,^h)f(y)dyd^ 

(2.23) b(x,£,h)= .Jf ,, / /W^^^^^^^^) 

x p(x, fa;) (fo. 

From (2.23) and a G S~°°, it is clear that 6 G S~°°. Using the Taylor expan- 
sion in h in (2.23) and a G S^°°, one gets easily b G S^°°. Thus T h $ (-h 2 A g ) G 
f^ 00 , and since Th is self-adjoint for the volume form dvh given by (1.3), 
one has also <J?o( — h 2 A g )Th G The proof of our lemma is complete □ 

Using the Taylor expansion a(x + hmv, £, h) = ^x a i x i C; M an d 

(mv) a e lt ^ mv = (d^/i) a e lt ^' mv , we get from (2.23) that the symbol 6 admits 
the usual asymptotic development, 

(2.24) b(x, e, fc) ^ V(7i/i) a ^<7(T ft )(a;,e, fc) ^a(x, £, 

L — ' a! s 

The following lemma will be crucial in our analysis. 

Lemma 3. Let <£> G Cg°([0,oo[), and A h = h~ 2 (T h - T d)h )<5> (-h 2 A g ) . 
Then Ah belongs to £~i°° ■ Its principal symbol, ao(Ah), satisfies near £ = 0, 

a (A h )(x,o = (^\^\l(r'M-r' d (o) 2 ) 

(2-25) 

where Ric{x) and S{x) are the Ricci tensor and the scalar curvature at x. 
Moreover, let U be a local chart, K a compact subset of U and (p G Cq°(U) 
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such that <p(x) = 1 in a neighborhood of K; let a(x,£,h) — ^(/i/i) fc cifc(a;, £) € 
S~j°° be such that in this local chart one has A^^p = Op(a) + Rh with 
smoothing. Then, for all k and all x £ K one has a,k(x,0) = 0. 

Proof. Let xq G M and let e±, . . . ,e d be an orthonormal basis of the 
tangent space T X0 M. For x = {x-y, . . . , xj) € M. d , we identify x with ^ Xjej € 
T XQ M. Let s i y exp xo (sx) be the geodesic curve starting at xo with speed x. 
Then, for r > small, the map <p XQ : x \— > exp (x) is a diffeomorphism of the 
Euclidean ball \x\ < r on an open neighborhood U of xq, and the coordinates 
xj in U are called geodesies coordinates centered at xq. In these coordinates, 
one has x$ = 0, and (gij(0)) = Id. Let R be the Riemann curvature tensor 

at x = and R(j,k)(l,m) = ( R (^> d§^ ) 9§T I dt~ ) ■ Then the Ricci tensor and 
the scalar curvature at x = are given by 

(2.26) Rici^— ,— J =#icj 1 * = 23%j)(i, fc ) ) 5 = Z] ffic w- 

Moreover, one has in these geodesic coordinates (see [15], page 474) 

(2.27) djgi :in (0) = 0, dj <9fc#z, m (0) = -\R(i,j){ m ,k) ~ \R{i,k){m,j) 
or, equivalently, 

(2.28) g id (x) = 5 id + \{R{x,ei)x\e 3 ) + 0(x*). 
Consequently, one has 

(2.29) ^det(g)(x) = 1 - \Ric(x, x) + 0(x 3 ). 

From this formula, parity arguments, and 2cdT' d (0) = — J^^yjdy, we get 

(2.30) \B(0,h)\ =h d c d ^l + ^^Sh 2 + 0(h 3 )y 

Moreover, in geodesic coordinates, one has k(0, u) = Id = m(0, w) and p(0, v) = 
y/det(g)(v), and thus from (2.19), (2.29), (2.30) and (1.4), we get 

<r(T h )(0,Z,h) = [ e*Vdet(s)(H<fr> 

\B{0,ti)\ 7| V |<1 

= T d (|e| 2 ) f 1 - - / e *«ffic(t>, «) d» 

V 3 J 6c d y|„|<! 

(2.31) +0(/i 3 ) 

= r j( l?l 2 ) + (-r,(|«|> )Sf 5 + i £ <«>) 

+ 0(/i 3 ). 
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Since G d (0 = T d (\£\ 2 ), one has 

-(0 = 2^, fc (r' d (o) + |ei 2 r'j(o)) + 4^e*r3(o) + o(|£| 4 ), 



9^- <9£ fc 

and from T d (\£\ 2 ) = 1 + r^(0)|£| 2 + 0(|£| 4 ), we get from (2.31), 
<r(T h )(0,£,h) 

(2.32) 

+ h 2 (^f(r d (o) - r' d (o) 2 ) + o + o(iei 4 )) 

+ 0(h i ). 

Let us now compute the symbol of the operator T d ,h^o{— h 2 A g ). Until the 
end of the proof we use the Einstein summation convention. First we remark 
that in local coordinates the symbol of the operator —h 2 A g is given by p = 
po + hpi with p {x,C) = g jk (x)CjCk = and pi(x,£) = -igkCk- Here (g jk ) 
denotes the inverse matrix of the matrix (gjk) and gu = d X] gi k + 2^Sp k d x .g 

where g is the determinant of the matrix (gjk)- Let F = &oT d and F be an 
almost analytic extension of F. Then 

(2.33) F(-h 2 A g ) = - [ dF(z)(-h 2 A g -z)- 1 L(dz), 

where L(dz) = dxdy is the Lebesgue measure on C and d = \ (d x + id y ). Let 
9? E Cq° be equal to 1 near x = 0. For any z£C\K there exist symbols 
ao,a\,a2 such that in local geodesic coordinates we have 



(2.34) (-h 2 A g - z) Op(a + hat + h 2 a 2 ) = <p(x) + h 3 R h 

'el' 



with Rh £ From the symbolic calculus it suffices to set 



Po- z Po ~ z J 

(2.35) _ i 

02 = (pottioi +Pofoa + piai +pittia ), 

Po- z 

where for two symbols /, g we define f^g(x, £) = E| Q |=j jh. d ff( x i d x*g( x , 
£). It follows that 

(2.36) F(-h 2 A g )ip(x) = Op(b + ^ + h 2 b 2 ) + h 3 R h 

with bj(x,£) = 7 J c dF(z)cij(z, x,£)L(dz) and G In particular we have 
bo = ip(x)F(\£\ 2 ), and as a±(z,0,^) = 0; we get also &i(0,£) = 0. Let us com- 
pute a 2 (z, 0,£). First, we observe that pi(0,£) = 0. Moreover, as d Xf ,po(0, £) = 
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0, for all k we have also (pi$ia )(z,0,g) = 0, poMi (2,0,0 = °( JlUlp) and 
PohaoMZ) = A9 (}$-$ m - Therefore, from (2.27) we get 

h(0,® = — / dF(z) 1 L(dz)Ag lm (0M m + 0(|g| 3 ) 

(2-37) = -^F"(|e| 2 )A<? im (0)^e m + 0(\Hf) 

= ^F"(0)RicU l U + O(\£\ 3 ). 

Therefore, we conclude that in geodesic coordinates, the symbol of F(—h 2 A g ) 
satisfies 

a(F(-h 2 A g ))(0,£,h) 

(2.38) 

= F(|£| 2 ) + h 2 H-P-Ricfa + 0(£ 3 )) + 0(h 3 ). 

Then, from (2.32), (2.38) and the rule of symbolic calculus, which are 
valid for Th by (2.24), we conclude that Ah belongs to £Z°° and that (2.25) 
holds true. 

Finally, since T h (l) = 1 = T d (-h 2 A g )(l) and <5> (-h 2 A g )(l) = $ (0), one 
has Ah(\) = 0; therefore Ah<p(x) = 0{h°°) for any x G K, and therefore, 
Op(a)(l)(x) = a(x,0, h) = 0(h°°) for any x G K. The proof of Lemma 3 is 
complete. □ 



The following lemma will be used in the sequel to handle the very high 
frequencies. 



Lemma 4. Let x G Co°W be equal to 1 near 0. There exists h > 0,C 
such that for all p G [1, oo], all h G ]0, ho] and all s > 1, one has 

(2-39) Th{ l- x )(zH^L\ 



Li' 



Proof. Set h = h/y/s. Then x(~~ R ) is a /i classical pseudo-differential 
operator, and belongs to the class £~ L . Let Rn(x, y) d g y be the kernel of the 
operator x(— 7— £ ). Then Rfi(x,y) is a smooth function of (x,y) G M x M, 
and for any a, there exists a nonincreasing function ^ Q with rapid decay 
such that for all H G ]0, 1], one has 



(2.40) \V% y R h {x,y)\<h- d -\ a ^ a 



d g (x,y) 
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Let Qh,s(x,y) d g y be the kernel of the operator T/j(l — x) ( ~ h ) • Then one 
has 

(2.41) Q h>s (x,y) = - ^^r^yr / Rn(z,y)d g z. 

\B(x,h)\ \B{x,h)\ J B(X)h) 

By the Shur lemma, it is sufficient to prove that there exists /io > 0, Cq such 
that 



xeM,he]o,h ] 

(2.42) 



sup / \@ htS (x,y)\d g y < C /y/s, 
\he]0M] J 

sup / \Q hjS (x,y)\ d g x < C /Vs. 
;he]o,hn\J 



yeM,he]0,ho] 

We shall prove the first line in (2.42), the proof of the second line being the 
same. One has h<h for s > 1, and from (2.41) and (2.40), we get that for 
any given Co > 0, one has for all h € ]0, co/2], 

d g (x,y) > c 

(2.43) 

=> \&h, s (^y)\ < K- d Mco/2h) e 0(tf°) c 0(s-°°). 

Thus we may work in a local chart U centered at a given xq E M, with local 
coordinates x = (x±, . . . , x^) E K rf , and we are reduced to prove in this local 
chart, for some Co > independent of xq, h E ]0, ho],s> 1, 

(2.44) sup / \Qh,.(xo = 0,y)\d g y<Co/^. 
he]0,ho]J\v\<2co 

Let/,(,) = 1 ^f> .Onehas 

(2.45) ©m( x o,!/) = f X0 (y) - / t R h (y,z)f Xo (z)d g z. 



Let r H (y,Z, h) ~ £ fe G be the symbol of % = x(-^ 2 A g ) E 

£ C 2°° in the local chart U. Then all the r^k{y,0 are smooth functions of 
(y,£) with support in < tq if is supported in r < tq. Moreover, by 
(2.11), one has rn,o(y, 0) = 1. Therefore, we get with bo(y,u) smooth in y 
and in the Schwartz class in u, and for some tp with rapid decay, 



(2.46) 



Rh(y, z)^detg(z) = H \ ( y, — ^— ) + q h (y, z) 



b (y,u)du = l, \q h (y,z)\<h- d+1 ip( l 
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Set y = hy, z = hz and 0^(0, y) = h d Oh tS (xo, y). Then (2.44) becomes 



(2.47) sup 

h£]0,h ]J\y\<'2coh-^ 

One has by (2.46), 



I ©h,s(0, y) I V det g (xq + hy) dy < C /y/s. 



(2.48) 



and 



(2.49) 



qh(y,z)f X0 (z)d g z 



< Ctr d+1 



d+lj, ( \y~ Z \ \ 1 {d 3 (x ,z)<h} 



J \B(x ,h)\ 



dz 



d g (0,hz)<h 



[ y^~Y~ ) 

y- z 
h 



h-% y 



(fx (y) - f xo (z))dz. 



From (2.48) and (2.49), we get for some with rapid decay, 
\Qh,.(P>y)\<cJ s d l 2 ^s\y-z\) 

x {ttt-{d g (Q,hz)<h} + \^{d g (Q,hz)<h} ~ 1 {d s (0,hy)<h}\)dz- 



(2.50) 



This implies 



L 



\y\<2c h- 



(2.51) 



<C 



<c h + 



|0h,s(O,y)| v / detc/(xo + hy)dy 



s d /^(^S\y-z\) 



x {ft*-{d g (0,hz)<h} + \^-{d g {0,hz)<h} ~ 1 {d g (0,hy)<h} I) dy 

—±du)<C /VS. 



1 POO 



'0 Jui/il + V 

The proof of our lemma is complete. □ 
3. The spectral theory of T^. 

3.1. Estimates on eigenf unctions. In this section, we prove estimates on 
the eigenfunctions of 2V Let us recall that ||/||_f/ s (M) denotes the usual 
Sobolev norm, and that the semi-classical Sobolev norm ||/||/i S is defined 
by (2.6). For a family f h E L 2 (M), we shall write f h £ Oc^ih 00 ) iff there 
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exists ho > 0, such that for any s,N there exists C S) n such that one has 
\\fh\\Hs(M) < Cs,Nh N for all h G ]0, ho]. If f h = £ fj.h^j is the Fourier expan- 
sion of fh in the basis of eigenfunctions of A g , this is equivalent to 

(3.1) ' 

Vj,V/iG]0,/i ]- 

Let < S < 1 and /to > 0. For /i G ]0, /io] , let e h be an eigenfunction of 
with ||e h || i 2 = 1, associated to an eigenvalue Zh G [6, 1], so (T^ — z/ l )e /l = 0. 

Lemma 5. There exists ho > 0, and for all j G N i/iere exists Cj > 0, 
such that, the following inequality holds true 

(3.2) sup \\e h \\ h j<Cj. 

he]o,h ] 

Proof. We use the notation of Lemma 2 and we set T^f. = ijj^Th^Pk- 
One has for h small enough = ^fc^/ifc- F° r am/ given k, we denote by 
x = (x\, . . . ,Xd) local coordinates in U k , and we choose a partition of unity 
in R d of the form 

(3.3) i=£*«, e a (x) = e(^^\ 

with 9 G Cg°. Then, for any integer m, there exists D m independent of h 



such that for any u G i/ m (M ) with compact support, one has 



(3.4) D" 1 ^ \\0au\\l m < \\u\\l tm <D m ^2 P 



i„\\ 2 

a a \\h,rn- 



If 9' G Cq° is equal to 1 on the set {X, dist(A A , support(#)) < 2}, one has 
for h g]0, ho] with ho > small enough, 9 a Th = 9 a Th9' a for all a. For any 
given a, we perform the change of variable x = h{a + X). Let S a be the 
rescaled operator acting on functions of the variable X defined by [with 
f(x)=F(^)] 

(3.5) 9 a T htk 9' a (f)(h(a + X)) = S a (F)(X). 

Let us first show that S a is the sum of two quantized canonical transforma- 
tions of degree —(1 + d)/2< —1. From the definition (3.5) of S a and (2.16), 
one has 

(3.6) x J e i( - x - Y ^q(h(a + X),Z,h)9'(Y)F(Y)dYd£, 
q(h(a + X),£,h)= [ e^ m ^ a+x ^ hv >p(h(a + X),hv)dv. 

J\v\<l 
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Let us compute the integral which defined q(x,£,h) for |£| large. The phase 
v —7- £.m(x, hv)v as no critical points in v, so if x{ r ) £ C^°]0, 2[ is equal to 1 
near r = 1, one has 

q(x, £, /i) = C x{r)r d ~ l ( f e^ rm{x ' hrLu)u) p({x, hru)) dJ) dr 
Jo \7|^|=i / 

(3.7) 

+ n(x,£,h), 

where n is a symbol in S~°°. The phase uj — > £.rm(x, hroj)u has two non- 
degenerate critical points on the sphere \u)\ = 1, oj^t = dh i^_i/ 2 Sl + 0(h), 
since ±7^—77-7-777 are the two nondegenerate critical points of the phase 

|s -1 / J (x)£| 

a; —7- £.rm(x, 0)w, and the critical values (homogeneous in £ of degree 1) are 
&±(x,r,£,h) = =b"|£| z + 0(/i) since |g _1//2 (x)<$;| = Using the stationary 
phase theorem, we get 



(x,£,h) = j\{r)r d -\e i ^ x ^h) a+{x ^ h) 



(3.8) 







where cr± are two symbols of degree — (d— l)/2. By integration in r, we thus 
get 

q(x, e, fc) = e »+(*.l,e,h) T+ ( Sj £ fc) + e »-(*,i,e,fc) r _ ( Sj £ h ) 

(3.9) 

+ n(x,£,/t), 

where r± are two symbols of degree — (d + l)/2. From (3.9) and (3.6), 
we get that S a is (uniformly in a,h for h £]0,ho] with ho > small), 
the sum of two quantized canonical transformations of degree — (d+ l)/2, 
with canonical relations closed to the ones associated to the phases (X — 
Y)£ ± |elh(a+X). that is, of the form (Y, 77) ^ (X = V ± ///M^ + = 

Since Th is (in the variable X) the sum of two quantized canonical trans- 
formations of degree —(1 + d)/2, and since e h = — T^e^), and z/i > <5, we 
get that there exists c and for all m, C m , independent of h,a, such that 

(3.10) \\0(X)e h (h(a + X))\\ Hm <C m \\e'(X)e h (h(a + X))\\ H . 



m— 1 . 

X 



where denotes the Sobolev space in variable X , as soon as 0'(X) is equal 
to 1 at each point X whose distance to support(#) is less than c. From (3.10) 
with m = 1, (3.4), and hd x = dx, we get for \{ x ) £ Co°(£4) and h e]0,/io] 
with h$ > small, 

(3.11) \\x(x)e h (x)\\ hA <C\\e h (x)\\ L 2 iUk) . 
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Therefore, since (£/&) is a covering of M, we get \\e \\h,i 

<C\\e h \\ L 2. We can 
now iterate this argument from (3.10), and we get for any j, 

(3.12) \\e h \\ h , 3 <C 3 \\e h \\ L 2. 
The proof of our lemma is complete. □ 

Remark that there exists s\ > 1 such that |r^(s)| < | for all s > si — 1. 
Let x S Cq°(M + ) be equal to 1 on [0, s±] and equal to on [s\ + l,oo[. 

Lemma 6. Let e h as in Lemma 5. Then 

(3.13) X {-h 2 ^ g )e h -e h = Oc^{h^). 

Proof. Let (ej)j g N be an Hilbertian basis of L 2 (M,d g x) such that 
—A g ej = Xjej and consider II S the orthogonal projector on spanje.,-, h 2 Xj > 
s}. By Lemma 4, there exist so,/io such that 

(3.14) Vs>s sup \\U s T h U s \\ L 2 < 5/2. 

he]0M] 

Let S2 > max(si + 1, so) and let X2, X3 De smooth functions such that 1r, = 
X + X2 + X3, X3(s) = for s < s 2 - 1 and x 3 (s) = 1 for s > s 2 . Let X2 G 
C Y o°(]R + ) equal to 1 near [s±, S2] and equal to on [0, s\ — 1] U [s2 + 1, 00 [. On 
supp(x2(s)) we have Zh — ^d(s) > f • Hence it follows from Lemma 3 that 
there exist E £ £° such that E(T h - z h ) = X2(~h 2 A g ) + R with R £ h 00 ^ 00 . 
As (T h - z h )e h = 0, we get 

(3.15) X2{-h 2 A g )e h £Oc°°{h°°). 
Set e h = Y, j x h j e j . Then 

u S2 e h -x3(-h 2 A g ) e h = j Y,u"i 2 ^M;<j 

h 2 Xj> s 2 j 

= - Yl Xsi^X^xfe. 

si<h 2 Xj<S2 



(3.16) 



As X2 = 1 on [si,S2], it follows from (3.15) and (3.1) that one has H S2 e h — 
X 3 (-h 2 A g )e h G C oo(h°°). Therefore we get 

(3.17) e h = x(-h 2 A g )e h + U S2 e h + O c ~ (h°°). 

Since II S2 is bounded by 1 on L 2 , applying H S2 (Th — Zh) to this equality, we 
get 

(3.18) U S2 (T h - z h )U S2 e h = -U S2 (T h - z h ) X (-h 2 A g )e h + O^h 00 ). 
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Let x £ Co°([0;Oo[) be supported in [0,S2[ and equal to 1 near the support 
of x- Then, thanks to Lemma 2, we have 

(3.19) (T h - z h ) X (-h 2 A g ) = x(-h 2 A g )(T h - z h ) X (-h 2 A g ) + h 00 ^ 00 . 

From (3.18), (3.19) and U S2 x(-h 2 A g ) = 0, we get U S2 (T h - z h )U S2 e h G 
Ojji{h°°). Since S2 > so, the operator H S2 (Th — 0/i)LT S2 is invertible on the 
space U S2 (L 2 (M)). Consequently, U S2 e h is 0(h°°) in L 2 (M). On the other 

hand, from Lemma 5, (3.2), one has for any integer j, \\Ag' 2 Il S2 e h \\^2 = 

\\H S2 A J g 2 e h \\i J 2 < Cjh~i . By interpolation it follows that for all j, one has 

\\A j / 2 Il S2 e h \\ L 2 G 0(h°°), that is, one has U S2 e h G O c ^{h°°). Then (3.13) 
follows from (3.17). The proof of our lemma is complete. □ 

For Zh G [5, 1], set Zh = l - h 2 Th, so that e h satisfies The h = (1 — h 2 Th)e h . 
The next lemma is a refinement of Lemma 5. 

Lemma 7. For aZZ j G N, t/iere exists Cj such that for all h G ]0, /to] j ^ e 
following inequality holds true: 

(3-20) \\e h \\ HHM) <C j (l + r h y/ 2 . 

PROOF. By Lemma 6, we have e h — x(—h 2 A g )e h G Oc°°(h°°), and there- 
fore using also Lemma 1, we get ((T h - l)x(-/i 2 A 5 ) + h 2 r h )e h G O c °°(h°°) 
and it follows from Lemma 3 and (T d - 1)(1 - x)(-/i 2 A 9 )e h G 0c->(ft°°) 
that 

(3.21) ((r d - l)(-h 2 A g ) + /i 2 ^ + h 2 r h )e h G O c =o (ft 00 ) 

with Ah G One has (T d — l)(s) = — sF^(s) with F d smooth, and from 

(3.21) , we get 

(3.22) - A g F d (-h 2 A g )e h = (A h + T h )e h + O c ™(h°°). 

Since is uniformly in h bounded on all H^{M), and ||e ft ||^2 = 1, we get 
from (3.22) for all j G N, with Cj independent of h, 

(3.23) \\F d (-h 2 A g )e h \\ Hj+ 2 {M) <Cj(l+T h )\\e h \\ 

Hi (M) ■ 

Since F d (s) on [0, s\ + 2], we get (3.20) by induction on j from (3.23) 
and (3.13). The proof of our lemma is complete. □ 

3.2. Proof of Theorem 1. Let us recall that there exists 70 < 1 such that 
T d (s) G [-7o,l] for all s£l Let e G ]0, (1 — 7o)/2[ and x (t) € Cg°([0, oo[) 
equal to 1 near £ = and such that x(t) G [0, 1] for all t. Thanks to Lemma 
4, there exists s > such that 



(3.24) 



T h (l - x) 



-h 2 A g 



L 2 (M,du h ) 
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On the other hand, thanks to Lemma 3 we can apply the Garding inequality 
to the pseudo-differential operator Thx( fc a Ag ) to get for h > small enough, 

(3.25) lT h x (^^) f, f) > (-70 - s)\\f\\ L 2 iM>dVh) , 

\ V s // L 2 {M,du h ) 

where we have used the fact that sup^ ||/||L 2 /ll/llL 2 (A/,d^) S oes to 1 when 
h goes to 0. Combining equations (3.24) and (3.25), we obtain 

(3.26) (T h f,f) L 2 {M4uh) > (-7o-2 £ )||/||| 2(MAfe) 

which proves the first statement of Theorem 1 as Ty t is self-adjoint on 
L 2 (M,du h ). 

Let us now prove (1.8). Set |A/J = 2(d + 2) 1 ~ f J h . For k < L, we denote by 
m k = dim(Ker(A g + X k )) the multiplicity of A^. Let po G Cg°(M) be equal 
to 1 near zero. Then there exists Hq > such that for h g]0, ho], one has 
e = po(—h 2 A g )e for any e G Ker(A g + A&) with k < L. Thus, if (£/j) is a 
finite covering of M by local charts and 1 = fj a partition of unity with 
(fj G C 3O (L r j), one has 

(3.27) (T h - T d , h )(e) = Y,{Th ~ T dth )p (-h 2 A g ) Vj (e). 

j 

From Lemma 3 one has for each j, (Th — T d ^h) P o(—h 2 A g )(pj = h 2 Op(a) + Rh, 
with a = ci2 + hd3 + • • • G S^ 00 compactly supported in x G Uj, Rh smooth- 
ing, ci2(x, £) = 0(£ 2 ) near £ = and 03(2;, 0) = 0. As e is smooth and does 
not depend on h, it follows that ((Th — T d ^h)po(—h 2 A g )ipj(e) G Ojji (h A ). 
Therefore, 

(3.28) \\(T h -T d , h )(e)\\ L2{MtdVh) = o(h 4 ). 

Moreover, F dA e = T d (h 2 X k )e = (1 + h 2 T' d (0)X k + 0(h A ))e. Combining this 
with (3.28) we obtain \\(\A h \ - h)e\\L 2 (M,dv h ) = 0(h 2 ) for all e G Ker(A g + 
Afc), and since |A^| is self-adjoint on L?(M,dvh), we get that there exists Co 
such that 

V/iG]0,/i ],V0<fc<L 

(3.29) 

card(S;pec(| A h |) n [X k - C h 2 , \ k + C h 2 ]) > m k . 

Now, if e h is a normalized eigenfunction of |A^|, |A/j|e ft = The h , with Th 
bounded, one has, by Lemma 6, e h — po(— h 2 A g )e h G Oc°°(h°°), and also 
by Lemma 7 since Th is bounded, ||e' l ||#j (M) < Cj for all j, with Cj inde- 
pendent of h. Thus the same argument as above shows that there exists C 
independent of h such that 

(3-30) \\{r h + ^ 9 )(e h )\\ L 2 {M4uh) <Ch 2 , 
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and thus dist(r/ l , Spec{— A g )) < Ch 2 . It remains to prove that for h small, we 
have equality in the right-hand side of (3.29). Let p > and let ei(h), . . . , e p (h) 
be a family of eigenfunctions of |A/J associated to the eigenvalues Tj(h) E 
[Afc — Coh 2 , Afc + Co/i 2 ], orthonormal for the scalar product (•, -)l 2 (m ,dv h )- By 
Lemma 7, there exists a sequence (h n ) going to zero as n — > oo such that 
£i{hn) converges in H 2 . Denoting fi its limit we get from (3.30), —A a fi = 
Xkfi for all I = l,...,p and the functions fi are orthogonal for the scalar 
product (■, •) L 2 (M,d g x)- This proves that to& > p, and completes the proof of 
(1.8). (In particular, this implies that 1 is a simple eigenvalue of T^.) 
Let us now prove the Weyl estimate (1.9). 

Let S e]0,l[ be given. Let r G [0,(1 - 5)h~ 2 ]. Observe that N(l-rh 2 ,h) 
is the number of eigenvalues of | A/J in the interval [0, 2{d + 2)r]. We denote 
by No(a,h) the number of eigenvalues of r^(— h 2 A g ) in the interval [a, 1]. 
Let us define the function &h(s) and the operator |A°| by the formulas, 

^(,) = 2(d + 2) 1 ~ r 2 d(s) , 

(3.31) 

\A° h \ = $ h (-h 2 A g ). 

Then Nq(1 — rh 2 ,h) is the number of eigenvalues of |A^| in the interval 
[0, 2{d + 2)t] . Let us first show that iVo satisfies the Weyl estimate (1.9), 
that is, there exists C such that for all h £ ]0, /io] an d all r € [0, (1 — 8)h~ 2 ], 
one has 

N (l-Th 2 ,h)-(2irh)~ d [ dxdi 

■/r d (|C|J)6[l-Tfta,l] 

< C ( 1+r) (<*-l)/2_ 

To prove this point, let n + (A) [resp. n _ (A)] be the number of eigenvalues 
Xj of — A g in the interval [0, A] (resp. [0,A[). By the classical Weyl estimate 
with accurate remainder (see [7]), one has 

(3.33) n ± (A) = (2vr)- d / dx d£ + 0(X^ d ~ 1 ^ 2 ). 

J\m<* 

By (3.31), Nq(1 — rh 2 , h) is the number of eigenvalues Xj of — A g such that 
1 - T d (h 2 Xj) < rh 2 . Since r < (1 - 5)/i~ 2 , the set {s > 0; 1 - r d (s) < r/i 2 } 
is a finite union of disjoint intervals Jq U ■ ■ ■ U 4 with Jo = [0, so(t/i 2 )], L- = 
[s~(r/i 2 ), sjirh 2 )] for 1 < j < k, and such that Co < s7/ < < < < 
• • • < sjjT < ci with Co > independent of /i, S and ci independent of h. Thus 
we get 

(3.34) N (l - rh 2 ,h) = n+(s h- 2 ) + ^n+(s+/i" 2 ) - n~(sjh~ 2 ). 

3=1 



(3.32) 
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Observe that k = when r < ch~ 2 with c small enough, and in that case 
one has by (1.6), soh~ 2 ~ 2(d + 2)r, and therefore (3.32) is consequence of 
(3.33). On the other hand, in the case r > chT 2 , then both (sfh' 2 )^- 1 ^ 2 

and (soh~ 2 )( d ~ 1 ^ 2 are of order t^ -1 )/ 2 , and thus we get (3.32) from (3.33) 
and (3.34). 

Let E T be the finite dimension space spanned by the eigenfunctions Bj 
of -A g with <f> h (h 2 \j) < 2r(d + 2). Then by (3.31), one has dim(.E T ) = 
N (l - rh 2 ,h). By (2.30) and |||A h ||| L 2 <Ch~ 2 , one has for all f£L 2 , 

(3.35) \(\A h \f\f) LHM!ZhdVh) - (\A h \f\f) L2(M4gX) \ < C"||/||| 2 . 

Let x S C^°([0, oo[) equal to 1 near the compact set {s > 0; 1 — T d (s) < 1 — 5}. 
Then / = x(~h 2 A g )f for all / G E T , and from Lemma 3, one has (jA/J — 
\ A °h\)x(-h 2 A g ) = -2(d+2)A h . Thus, since A h is bounded on L 2 , from (3.35) 
we get that there exists C_ = C_(<5) independent of r,h, such that for all 
f € E T , one has 

(3.36) {\A h \f\f) L , {NLZhdVh) < 2(r + C_)(d + 2)||/||| Whdl/h) , 
and this implies, by the min-max, 

(3.37) N (l - rh 2 , h) = dim(£ r ) < N(l - (r + C-)h 2 ,h). 

Let .F T be the orthogonal complement of £7 T in L 2 (M,d g x). Let G Cg° 
such that ||7X(1 - 9){-h 2 A g )\\ L 2 < 5. Let x G Cq with values in [0, 1], equal 
to 1 near [0, 1 — 6] Usupport(#). Let tp = 1 — x, so that (1 — Q)ip = tp. Let = 
(\A h \ - \A° h \) X (-h 2 A g ) G «5 C 7°° and B h = X (-h 2 A g )(\A h \ - |A°|) G f J 00 be 
given by Lemma 3. Then, one has 

(3.38) \A h \ = x\A° h \ x + i>\Ah\x + xK|^ + ^A^ + A h + B h ^. 

The operator Ah + is bounded on L 2 by a constant C(S) uniformly in 
h. Prom tp(l — Th)ip = ip 2 — tpT^l — 9)ip, we get 

(3.39) {MA h \i>f\f) LWgX) > 2(1 -5)^Uf\\h {M , da *y 

Therefore, from (3.35) we get that there exists C + = C+(5) > independent 
of t, h, such that for all / = £^ >T x j e j G -F T , one has 

(\A h \f\f) L2(M)ZhdVh) + (d + 2)C + \\f\\ 

> Y^^X^ix 2 + ^){h 2 \ J )\x J \ 2 

\j>T 

(3.40) 

Xj>T 

>2r(d + 2) ^ k i | 2 >(2r(d + 2)-C/ l 2 )||/||i 2(MiZ;i(il/h) , 
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and this implies by the min-max for r large enough, and h g]0, ho] with ho 
small, 

(3.41) N (l - rh 2 ,h) = codim(F r ) > N(l - (r - C+)h 2 , h). 

Then we obtain the Weyl estimate (1.9) from (3.32), (3.37) and (3.41). Fi- 
nally, the estimate (1.11) is an easy byproduct of the estimates (3.20) of 
Lemma 7. The proof of Theorem 1 is complete. 

3.3. Proof of Theorem 2. Let us recall that &h(s) and |A°| are defined 
in (3.31). 

One has 2(d + 2)(1 — T^s)) > c\ min(s, 1) with c\ > 0, and, therefore, 

(3.42) $h(h 2 \j) > ci min(Aj, h~ 2 ). 

Observe that there exists ho,co > such that for all z Gil, all h £]0,ho], 
and all j 6 N, one has 

(3.43) \z-$ h (h 2 \j)\ >co(l + |«| +min(Aj,/i- 2 )). 

To see this fact, just observe that by (3.42), for c± min(A./, h~ 2 ) > A+l, 
(3.43) holds true, since z EU. Now, c\ mm(Xj , h~ 2 ) < A + l implies if 
is small, Xj < (A + l)/c\, and therefore, \$h(h 2 \j) — \j\ < C2h 2 , and (3.43) 
holds true also in that case since z £U. Since for h 2 Xj < C3 with C3 > small, 
one has \&h{h 2 \j) — Xj\ < c^h Xj, we get from (3.43), that there exists C 
such that for all z E U and all h 6 ]0, ho], one has 

1 1 



(3.44) sup 



<C7i 2 , 



z - <& h (h 2 Xj) z-Xj 
and this implies, obviously, 

(3.45) ll^-IA^-^ + A,)- 1 !^^^ 2 , 
and thus we are reduced to prove the estimate 

(3.46) ||(z - lA^I)- 1 - - lAOl)" 1 !!^ < C^ 2 . 

Observe that, as a straightforward consequence of Theorem 1 and of the self- 
adjointness of |A^| and |A°|, respectively, on L 2 (M,duh) and L 2 (M,d g x), 
there exists C > and ho > such that for all z £ U and all h E ]0, ho], 

(3.47) IKz-IA.mb + IKz-KiniL^ 



1+2 



Therefore, in order to prove (3.46), we may, and will assume that z satisfies 
h 2 \z\ < a, with a > small. Using Lemma 4, we then choose xo £ Co° ec L ua l 
to 1 on [0, so]) with support in [0, 2so], and, such that, 

(3.48) \\2(d + 2)T h (l- X o)(-h 2 A g )\\ L 2<d + 2-a/2. 
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Let x G Co° equal to 1 near [0,3s ], and set R h = (z- \A\ h )~ l - (z- |A|£) _1 . 
Then since |A^| commutes with A g , one has 

R hX (-h 2 A g ) 

(3.49) 

= (z- \A h \)-\\A h \ - \A° h \)x(-h 2 A g )(z - \A h \r\ 

From Lemma 3, one has (\A h \ - \A° h \)x(-h 2 A g ) = A h x'(-h 2 A g ), with x' 
equal to 1 near the support of %■> an d the operator E £~i°° satisfies 

(3-50) \\A h f\\ L 2 {M) <Ch 2 \\f\\ HHM) . 

On the other hand, from (3.43), we get 

(3.51) \\ X '(-h 2 A g )(z - \A° h \)- l f\\ H 2 iM) < C||/|| i2(M) . 
From (3.47), (3.49), (3.50) and (3.51), we get 

(3.52) \\R hX (-h 2 A g )\\ L 2<Ch 2 . 

It remains to estimate i?h(l — x)(~ h 2 A g ), and it is obviously sufficient to 
prove the two estimates 

(3.53) ||(* - lA^ir^l - X )(-h 2 A g )\\ L 2 < Ch 2 , 

(3.54) ||(* - KD^l - X )(-h 2 A g )\\ L2 < Ch 2 . 

Since x( s ) = 1 near s = 0, (3.54) is a consequence of (3.43). Let g G L 2 {M) 
with \\g\\ L 2 = 1 and let f = (z - lA^-^l - X )(-h 2 A g )g. Then * 

(3.55) (h 2 z - 2(d + 2)(1 - T h ))f = h 2 (l- X )(-h 2 A g )g. 

Let xi> X2 G Co° with support in [0, 3so[, with xi equal to 1 near the support 
of xx- One has Xi(l — x) = 0, and thus, multiplying (3.55) by Xi(~ h 2 A g ) 
and using Lemma 3, we obtain 

(3.56) h\z - \A\\)x\{-h 2 A g )J = h 2 A h x 2 {-h 2 A g )f + O c ™(h°°). 

Since on the support of xi> one has h 2 Xj < 3so, we get from (3.43), (3.47) 
and (3.56) that one has ||xi(— h 2 A g )f\\ H 2 < C; thus, since xi 1S arbitrary, 
||X2 ( — h 2 A g )f\\ H 2 < C, and from (3.56) and (3.50), we thus get 

(3-57) \\ X i(-h 2 A g )f\\ L 2<Ch 2 . 

Then, we deduce from (3.55) and (3.57) 

(3.58) (h 2 z - 2(d + 2) + 2(d + 2)T ft (l - X o(-h 2 A g )))f G L 2 (h 2 ), 

and from (3.48), we get ||/||i2 < Ch 2 . The proof of Theorem 2 is com- 
plete. 
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4. Proof of Theorem 3. 

4.1. The spectral theory of the Metropolis kernel. In this section, we will 
deduce from the results of Section 3, useful properties on the spectral theory 
of the Metropolis operator Mh- Let us write 

(4.1) M h = T h + R h . 

Then from (1.16) and (1.17), one has 

Rh{f){x) = m h (x)f{x) 



(4.2) 



Jd g (x, y )<h \\B(y,h\ \B(x,h\ 



Let a(x,y,h) < be the function 

(4.3) a(x,y,h) = h d ~ 2 min( — — — - — — - ,0 



B(y,h\ \B(x,h\ 

Then a is a Lipschitz function in x and y, and from (2.30), we get that there 
exists C independent of x,y,h such that 

(4.4) \a(x,y,h)\<Cd g (x,y), \V x a(x,y,h)\ + \V y a(x,y,h)\ < C. 

Since i?/ l (l) =0, one has mh{x) = —h 2 ~ d f d ^ y ^ <h a(x,y,h) d g y, and there- 
fore the function is Lipschitz and satisfies ||77i/i||l«> < C/i 3 and || Vmh\\L°° < 
Ch 2 . From these facts, one easily gets that there exists C independent of 
p G [l,oo] and h such that 

\\Rh\\LP<Ch 3 , 

(4.5) 

\\Rh\\w^<Ch 2 , 

where W l,p = {/ G L p , V/ G L p } is the usual Sobolev space. Therefore, 
is a small perturbation of T^. In particular, there still exist ho > and 7 < 1 
such that the spectrum of Mh is a subset of [—7, 1], 1 is a simple eigenvalue 
of Mh and since Hm^Hioo < Ch? and m/ l (x) > 0, the spectrum of Mh is 
discrete outside [0,Ch 3 ]. Let 

(4.6) Ch 3 < ■ ■ ■ < Ji k+ i{h) <Mh)<---< Mh) < Mh) = 1 
be the decreasing sequence of positive eigenvalues of Mh ■ Set 

h 2 

(4.7) 1-M h = — r|AJ. 

v ; 2(d + 2)' 1 

Then from (4.5), one has 

(4-8) \\\Ah\-\A h \\\ L 2<Ch. 
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From Theorem 1 and (4.8) we get that for any given L > 0, there exists C 
such that for all h G ]0, ho] and all k < L, one has 

l-jl k (h) X k 



Moreover, since ||T/j — M/ l || i 2 < Ch 3 , the Weyl estimate (1.9) remains valid 
for the number N(a, h) of eigenvalues ofM/jinthe interval [a, 1]: for o~g]0,1[, 
one has 

iV(l - rh 2 , h) - (2Trh)- d [ dx d£ 

Jr d Qt\Z)e[i-Th a ,i] 

^^(l+r)^ 1 )/ 2 

for any t G [0, (1 — 8)h~ 2 ], and therefore, the estimate (1.10) is still valid; 
for any r G [0, (1 — 5)h~ 2 ], one has 

(4.11) N(l -rh 2 ,h) <Cs{l + T) d ' 2 . 

The main result of this section is to prove that there exist Cs such that 
for any eigenfunction ejj of Mh associated to the eigenvalue Jlk(h) G [8,1], 
the inequality (1.11) still holds true, that is, with Tk{h) = h~ 2 (l — Jlk(h)), 
one has 

(4.12) ||e£|U» <C s (l + T k (h)f 4 \\e£\\ L 2. 

We will obtain this estimate as a consequence of (4.5), using Sobolev in- 
equalities and the following lemma. 

Lemma 8. Let N > 1, p G [l,oo] and 8 G]0,1[. Let s > suc/i that 
\l"d( s )\ — 8/2 for s > sq. Let xo G Co° suc ^ Xo( s ) = 1 on [0, so]- There 
exist C,Cn, ho, and for all z G K = {z G C, \z\ G [5, 2]} and a// /i G ]0, ho], 
operators E z ^,M z ,h which satisfy 

(4.13) E Zth (T h -z) = l- Xo(-h 2 A g ) +M z , h , 
and such that the following estimates holds true: 

\\E z h\\ijp < C, \\E z h\\ w i, P <C, 

(4.14) 

\Wz,h\\LP <Cfth , \Wz,h\\w^P ^^Nh . 

Proof. Let x G Co°([0,2[) ec l ual to 1 on I ' !]> and set = x(t/ s )- 
By Lemma 4, there exist so such that for all s > sq, one has ||7) t (l — 
Xs(—h 2 A g ))\\iP < 8/2. We then take s > so such that Xs = 1 near the sup- 
port of xo; and we set tp = 1 — Xs and tp' = 1 — Xis- For z £ K, T^ip — z is 
then invertible on L p . Set 

(4.15) E 1 = ij'(T h ^-zr\ 
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Then, there exists C, ho such that for all h G ]0, ho] and all z G K one has 

(4.16) H^illiP + 11^111^1* <C. 

The LP bound is obvious since operators in £~i°° are bounded on LP and ip' = 
1 — Xis] let us prove the W 1,p bound in (4.16). We denote by B any operator 
which is, uniformly in h > small, and z G K, bounded on LP. Let X be a 
vector field on M. Then by (2.16), one has [T h , X] = hB x X + B 2 . Thus, with 
L = T h ip - z, we get [L,X] = hB 3 X + B A and [X,^ 1 ] = hB 5 XL _1 + B 6 . 
Since for h small, 1 — hB§ is invertible on LP , we obtain XLr 1 = B 7 X + B 8 , 
and thus (4.16) holds true, since E\ = ij/L . Let <p G Cg°([0,3s[); from 
ip'(/> = 0, we get E\L(p = 0, and therefore 

(4.17) E l (f) = E 1 [cp,L\L- 1 . 

By Lemma 3, one has [(f), L] G h£^°° . Thus (4.17) implies ||£'i0|| L p + 
< Ch, and since <\> is arbitrary, by an easy induction from (4.17), 
we get \\E x 4>\ L v + \\Ei(j)\\ w i, P < C N h N for all N. Thus one has 

(4.18) E 1 (T h -z)=^' +M 

withM =^T fc (l-V) = E 1 {(j)T hXs + O(h°°£~ co )) if = 1 near [0,2s]. Thus 
A/i satisfies for all JV, 

(4.19) ||M||iP + ||M|| W i*<CW^- 

Now, by the symbolic calculus, there exist E 2 G and A/2 G h ca £~^° such 
that 

(4.20) ^ 2 (T/ l -z) = X4,-Ao+AA 2 . 

Here we use Lemma 4 and the fact that Th — z is elliptic near the support 
of X4s — Xo- Then = E\ + E 2 and A4,/i = M\ + M 2 satisfies (4.13) and 
(4.14). The proof of our lemma is complete. □ 

Let us now achieve the proof of (4.12). Let Jl(h) G [6,1] and He^H^a = 1. 
Then (Mh — ^(h))^ = is equivalent to (Th — Jl(h) + Rh)e h = 0, and using 
Lemma 8, we get 

(4.21) (1 - Xo )? h + (M m , h + E m , h Rh)^ = 0. 

Set ei = Xo^ h ) and e+ = (1 — Xo)^) so that e h = ej + e+. Since by (4.5) 
and (4.13) the operator M^ h ),h + E p(h),hRh is 0(h 2 ) on L p and W 1 ^ , we 
can solve equation (4.21) for e + on the form 

e+ = Sjir h \ h (ei), 

(4.22) 

ll<%(h),hllLP + l|'%(h) ) h|liv 1 .i' — C^ 2 - 
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Let l-h 2 r = Ji(h) and u = y/T+T. One has \A h \(e h ) = 2(d+2)(r + h~ 2 R h )(e h ), 
and therefore, with (\Ah\ — |A°|)xo = Ah, we get the equation 

(4.23) | A° |xo(^) = (2(d + 2) X0 (r + /i" 2 ^) - A + [| A h |, Xo])(^)- 

By (4.5) and Lemma 3, the operator 2(d + 2)xo( r + h~ 2 Rh) — Ah + [| A^|, xo]) 
is bounded by O(co 2 ) on L p , uniformly in h. Then by (4.22) and (4.23), we 
get for some p± £ ]d,oo[ and all p G [2,p*], that the following estimates holds 
true, with C independent of h: 

\\e h \\Lv<Cu d / 2 - d lv, 

(4.24) 

\\^\\ w ,, v <Cu d ^- d ^ +1 . 

Indeed, by (3.31) and (3.42), for xi £ Co° ec l ua l to 1 near the support of xo> 
one has |A°|xi = —A g Bh with Bh £ £~i°° elliptic near the support of xo- 
Thus, 1 1 e^ - 1 1 = 1 and (4.23) implies ||ei||^2,2 < Cui 2 , and thus ||ei||^i,2 < 
Cuj, so using (4.22), one gets that (4.24) holds true for p = 2. This also shows 
easily that (4.24) holds true for d = 2. When d > 3, then if (4.24) holds true 
for some p £ [2,d[, then let q G]p,oo[ be defined by d/q = d/p — 1. Then 
the injection W 1,p C L q shows that the first line of (4.24) holds true for q. 
Moreover, in (4.23), classical properties of — A g and the fact that operators 
in £^°° are bounded on W s > p , shows that ||ei|| W 2, P < Coj d / 2 ~ d / p+2 . Then 
the injection W 2 ' p C W l,q and (4.22) implies that the second line of (4.24) 
holds true for q. Then, from (4.24), we conclude the proof of (4.12) by the 
interpolation inequality for p* > d, 

(4-25) IMIl~<C||u||*;5 p 1u||^- 

4.2. The total variation estimate. In this section, we prove Theorem 3. 
Let IIo be the orthogonal projector in L 2 (M,dfiM) on the space of constant 
functions 

Then 

(4.27) 2 sup \\M%{x,dy) - duMhv = \\M% - T1 \\ L ^ L ~. 

Thus, we have to prove that there exist A, ho, such that for any n and any 
h £ ]0, ho], one has 

(4.28) e -Y(h)nh* < || M n _ n || Loo ^ L oo < Ae-^ h)nh2 

with j(h) ~ Oh) ~ 2(fe) when h -> °- Since ( M h - n o)(e?) = (1 - h 2 T^) n e\, 
with \ti — 2 (d+2) 1 — ^ ^ (4-9), the lower bound in (4.28) is obvious, and 



SEMI-CLASSICAL RANDOM WALK 29 

to prove the upper bound, we may assume n > C§h~ 2 . Let 5 g]0, 1[ be such 
that the spectrum of is contained in [—5, 1] . Then write — Hq = 
M h>1 + M ht2 with 

M h , 1 (x,y)= (l-h 2 Tk(h))e>;(x)ei;(y), 

(4.29) 

M h>2 = M h -U -M hil . 

Here 1 - h 2 r k (h) = Ji k {h). One has M% - n = M% x + M£ 2 , and we will get 
the upper bound in (4.28) for each of the 2 terms. Prom (4.29) and (4.12), 
there exist some a > such that 

(4.30) \\MH,x\\l<*>^l<*>< E (l-h 2 T k (h)) n (l + T k (h)) a . 

Tl(h)<T k (h)<(l-S)h~ 2 

Using 1 — x < e~ x , and the estimate (4.11) on the number of eigenvalues of 
Mh in [1 — h 2 r, 1], one gets for some C, /3, 

/•oo 

(4.31) ||A#i||l<»_> l <»<C / e- nh2x (l+xfdx, 

Jn(h) 

and we get for some C", 

(4.32) HM^IIloc^oo <C'e- nh2 ^ Vn>C /i" 2 . 

Since M™ is bounded by 1 on L°° , we get from M™ - n = + M£ 2 
and (4.31) that there exist Ci,m such that ||M^ 2 ||£°o_;.£oo < C\h~ m for all 
n > 1. Next we use (1.15) to write = + JCh with 

< 7 < 1, 

(4.33) 

\\ic h \\ L ^ LOO <c 2 h- d / 2 . 

From this, we deduce that for any p = 1,2,... one has = A p ^ + B Pt h, with 
Ai t h = rrih, Bi h = K-h and the recurrence relation A p+ i^ = m^Ap^, -B p +i,/i = 
TnhB P) h + K-h,M^. Thus one gets since is bounded by 1 on L 2 , 

II^/iIIl^^l 00 < 7 P , 

(4.34) 

llBp.fcllia^cc < C 2 h- d ' 2 (l + 7 + • • • + 7 P ) < C 2 ^ d / 2 /(l " 7)- 

Observe that \\M£ 2 \\ L oo^ L 2 < \\M£ 2 \\ L 2^ L 2 < 5 n and for g,p> 1, one gets, 
using (4.34), 

(4.35) < \\A Pt f l M^ 2 \\L°°->L< x > + ll-Bp^-^^lli 00 ^^ 00 

< + C 2 h- d ' 2 5 q /{\ - 7), 
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and this implies for some C, fx > 0, 

(4.36) ||Af£ 2 || L «,_> £ <» < Ce" n ^ Vn > l/h, 

and thus the contribution of MJ^ 2 is far smaller than the bound we have to 
prove in (4.28). The proof of Theorem 3 is complete. 

APPENDIX: CONVERGENCE TO THE BROWNIAN MOTION 

The purpose of this appendix is to answer a question of one of the referees 
about the convergence of the previous Metropolis chain to the Brownian 
motion on a Riemannian manifold (M,g). One classical and efficient way 
to prove such convergence is the use of Dirichlet forms (see [9]). Here, we 
present a self-contained proof, in the spirit of ([12], Chapter 2.4), making 
use of our previous results. The two main estimates are: the large deviation 
estimate (A. 15) of Proposition 1, and the "central limit" theorem (A. 46) of 
Proposition 2. 

We refer to [8] and [11] for a construction of the Brownian motion on 
(M, g) . For a given xq € M, let X X0 = {bo E C°([0,oo[, M),w(0) = x } be the 
set of continuous paths from [0,oo[ to M, starting at xq, equipped with 
the topology of uniform convergence on compact subsets of [0, oo[, and let 
B be the Borel <r-field generated by the open sets in X XQ . Let W XQ be the 
Wiener measure on X Xo , and let pt(x , y) d g y be the heat kernel, that is, the 
kernel of the self-adjoint operator e*^ 9 / 2 . Then W Xo is the unique probability 
on (X Xo ,B), such that for any < t\ < t 2 < ■ ■ ■ < tf. and any Borel sets 
A\, . . . ,Ak in M, one has 

Wxoiufa) G Ai,u)(ti) eA 2 ,.. .,u(tk) G A k) 
(A.l) = / pt k -t k _ 1 (x k ,x k - 1 )---p t;2 -t 1 (x2,x 1 ) 

J AixA2X---xA k 

x p h (x\, x Q ) dgXi d g x 2 ■ ■ ■ d g x k . 
For h e]0, 1], let M^x ^ e the closed subset of the product space M N , 
(A.2) Mh lXo = {x = (x 1 ,x 2 ,...,x n ,...),\/j > 0,d g (xj,xj+x) < h}. 

Equipped with the product topology, M N is a compact metrisable space, 
and the Metropolis chain starting at xq defines a probability V XQ ,h ° n M^, 
such that 'Pxchi-M^ XQ ) = 1, by setting for all k and all Borel sets Ai,...,A k 
in M, 

V xo h{xx gAi,i 2 € A 2 ,...,x k eA k ) 

(A.3) 

= / M h (x k -i,dx k ) ■ ■ ■ M h (xi,dx 2 )M h (x ,dxi), 

J AixA 2 x---xA k 
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where the Metropolis kernel Mh(x,dy) is defined in (1.15). Let j Xo ,h 

be the 

map from into X XQ defined by 

(A.4) j xo ,h{x) = uo Vj>0 u(jh 2 /(d + 2)) = Xj 

and 

jh 2 (j + l)/i 2 " 



(A.5) 



uj(t) is the geodesic curve connecting 



d + 2' d + 2 
j to Xj + \ at speed h~ 2 (d + 2)d g (xj, Xj+i). 



Observe that for /i > given, smaller than the injectivity radius of the Rie- 
mannian manifold M, the map j XQ ,h is well defined and continuous. Let P Xo ,h 
be the probability on X XQ defined as the image of V xo h by the continuous 
map j xo ,h- Our aim is to prove that P xo ,h converges weakly to the Wiener 
measure W XQ when h — > 0. 

Theorem 4. For any bounded continuous function uj \— > f(uj) on X xo , 
one has 

(A.6) lim / fdP XOih = [ fdW X0 . 

Observe that the proof below shows that our study of the Metropolis chain 
on the manifold M is also a way to prove the existence of the Brownian 
motion on M. 

Let us recall that the Metropolis operator Mh acting on L? = L 2 (M, d\iM) 
with d/iM = d g x/ Vol(M)) is defined by (1.17). If (p is a Lipschitz function 
on M, we denote by Mft i¥ , the bounded operator on L? defined by 

M h)V = e^ h M h e~^ h . 

The first ingredient we use in the proof of Theorem 4 is the following lemma, 
which gives an L 2 -estimate on the resolvent (z — M^) -1 near z = 1. 

Lemma 9. Let ip be a real valued Lipschitz function on M , p > and 
< 9 < 2-7T. Let us assume that the following inequality holds true: 

00 p k/2 H\\T- 

(A.7) psm{9/2) - ^ fc 7 I sin((fc " l )°l 2 )\ = c > °- 

k=2 

Then, with w = pe ld G C \ [0,oo[ and ip = ip l l 2 e te l 2 il), one has 
(A.8) ||(1 - M h ^-w)~ l \\ L 2 <l/c. 
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Proof. If k(x,y) is a complex valued bounded measurable function on 
M x M, let Ak,h be the bounded operator on L 2 , 

(A.9) A Kh {f)(x)= [ mhl (—l ) —LJ\k(x,y)f(y)d g y. 
Jd g (x, y )<h \\B(x,h)\ \B(y,h)\J 

With k*(x,y) = k(y,x), the adjoint on L 2 of Ak,h 1S equal to Ak*,h, and one 
has the obvious estimate 

(A. 10) 1 1 -4*^1 1 £2 < \\k\\L°°(MxM)- 

From (1.15) and (1.2), one has Mh = rrih + Ai.h, and an easy calculation 
gives 

(All) M h)V = m h + A^, h , k v {x,y) = l dg(x , y) < h e Mx) -^ /h . 

Let r(x,y) = l dg {x,y)<hi{^{x) - ip(y))/h. With tp = ip 1 / 2 ^ 9 / 2 ^, we thus get 

(A.12) M http = m h + 22 fc , A",h- 

k=0 

From (A.12) and w = pe ld , we get with S = —e~ %e l 2 (\ — — w), 
S = - e - i6 / 2 (l - M h ) + p 1/2 A Th + pe ie ' 2 Id + N, 

(A.13) 

«/2^ (P 1/2 e" /2 ) * 
k\ 



k=2 

Since r* = r, the second term in the first line of (A.13) is self-adjoint, and 
we get 

Im(S) = sin(0/2)(l - M h ) + psm(6/2)Id + lm(N), 

(A.14) 

Im(iV) = EV S[n({k ~ WMMh- 

k=2 

From sin(0/2)(l — Mh) > 0, and since from (A. 10) the self-adjoint operator 
A T k^ h has norm < 1 1 1 1 Lips ' we S e ^ f rom (A. 7) and (A.14) that Im(5) > eld. 
The proof of Lemma 9 is complete. □ 



From Lemma 9, we shall now deduce a key estimate on the probability 
that X% Xo , the nth step of the Metropolis chain starting at xq, satisfies 
d g (X£ xo , xq) > e. Let £o > be smaller than the injectivity radius of the 
Riemannian manifold M. 
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Proposition 1. There exist positive constants C,A,a,CQ,ho>0 such 
that for all e e]0,£o]> a ^ $ £]0,co£ 2 ] and all h E]0, ho], the following in- 
equality holds true: 

(A.15) sup V Xo , h (d g (X^ X0 ,x ) >e)< Ce^eT^I 6 . 

x £M,nh 2 <8 

Proof. We may assume nh > e, since otherwise V XOi h(,dg(X^ x ,xq) > 
e) = 0. In the proof, we denote by a,A,C positive constants, changing from 
line to line, but which are independent of h, e, xq E M and n > 1. One has 



r X0>h (d g (X% X0 ,x ) >e)= M%(x ,dy) 

Jdg{y,x )>e 

(A.16) 

= M%(l dg{ytXo)>E )(x ). 

Let ip(r) E C°°([0, oo[) be a nondecreasing function equal to for r < 3/4 
and equal to 1 for r > 1. For e E ]0,£o] an d £ ^\ set 

(A.17) 9xo 4x) = ^(^p y 

Then v^ce is a smooth function, and from ld g (y,x )>e — V^o.e < 1; we get, 
since is Markovian, 

(A.18) M£(l ds(2/ja:o)>£ ) < M£(^ 0>E ) < M£(l) = 1. 

We first deduce from Lemma 9 the following estimates on MV"{tp xoe ). 

Lemma 10. There exists cq > such that for nh 2 < cqE 2 , the following 
inequalities hold true: 

(A.19) \\MJi( Vxo , s )\\ L2{BM2)) < Ce-°* V"' 12 ; 

(A.20) ||^(^o, e )IU-(B(x 0|e /4)) < Chr^e-^l^. 

Proof. By the Cauchy-Schwarz formula, the self-adjoint operator MJf 
is equal to 

(A.21) M% = ^jz n {z-M h y l dz, 

where a is a contour in the complex plane surrounding the spectrum of 
with the counter-clockwise orientation. Let 9q E]0,7r/2[ close to ir/2 and 
po > small be given. Since we know that the spectrum of Mh is a subset 
of [—7, 1] with 7 E [0, 1[, we may choose a in the form o\ U o~2, with 

a 1 = {z = l- w(e),w(0) = p(0)e ie ,O E [0 , 2vr - O }}, 
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where the function p(9) > takes small values, is such that p(9) = p(2ir - 
Po = p(@o) an d will be chosen later, and with q = |1 — poe %e °\ < 1, 

o- 2 C {\z\ < q,dist(z, [-7, 1]) > p o sin(0o)}. 



Il.vllr-' 



Set g = ip XOt£ and f z = (z- M h ) l g. For z G <r 2 , one has ||/ Z || L 2 < posin( g o) , 
and from [z - m h )f z = M,hfz + 9, \z - m h (x)\ > dist(z,[0,l]) > p sm(6 ), 
and ||^4.i,h/a||i°° < C/i _rf / 2 ||/^|| i 2 , we get for z G o"2, with a constant C 
changing from line to line, 



11/ 



s\\L° 



< 



1 



{\\M,hfz\\L°° + |bl|L° 



(A.22) 
This gives 

(A.23) 



Posin(6> 

<Ch- d / 2 ( PoS m(e Q )r 2 



z n {z-M h y\g)dz 
z n {z-M h )-\g)dz 



<Cq n (p sm(9 )r\ 



<Cq n h- d / 2 (p sm(e ))- 2 . 



Observe that since nh > e, one has q n < e~ a£ / h < e ~ a£2 / nfe2 . Thus (A.23) 
gives (A. 19) and (A. 20) for the contribution of o<i- Next we use Lemma 9 to 
bound the contribution of o\ in (A. 21). 

Let p < 1 and set ip(x) = pV2dist(x,B(xo,e/2)). One has 1 1 "0 1 1 Lips = 
pV2, and if p(9) > is small enough, inequality (A. 7) is fulfilled with 
a constant c~ p{6) sin(0/2)(l - p) + 0(p 3/2 (9)) ~ p{9). From (A.8), and 
(z - M hjip )e^ h f z = e^^g, we get for z = 1 - w(9) G ax, since 99 = on 
B(x ,e/2), g = on B(x ,3e/4), and |e^/ h | = | e ™ 1/2 (eW^dist(^B(x ,e/2))A| 5 



fz||i 2 (B(a;o,e/2)) < II e^/z || L 2 < 



(A.24) 



l^lb 



< 



—I 



Jw 



One has (z — mh)f z = Ai^ifz) + 9 with g = on B(xo,e/2), and /i < cqE 
since /ie < n/1 2 < cqE 2 . For cq < 1/4, we thus get from (A.24), 



(A.25) 



I z\\L°° 



(B(x ,e/4)) < Cp 1 (9)h 



-d/2\ iw 1 / 2 ae/hi 



On a±, we set z = 1 — w = 1 — u , u = p ' (9)e 



l/2 ( mJ9/2 



w 



l l 2 . Then one has 



(A.26) 



z n {z-M h )-\g)dz= / {l- ^Yh^udu 
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where 7 is a contour in the upper half plane Im(u) > connecting u- 



V2 -i9 /2 



-Po 



(A.27) 



to u + = pl /2 e ie °/ 2 . From (A. 24), (A.25) and (A.26), we deduce 



z n (z-M h )- 1 (g)dz 



z n {z-M h )-\g)dz 



< CJ, 



L*{B(x ,e/2)) 



< Ch- d/2 J, 



L°°(B{x ,e/4)) 

where J is defined by (with a > small) 

2\n Juae/h 1 I^H 



(A.28) 



|(1-^)V 



and it remains to verify that J satisfies 



(A.29) 



J < C 2 e~ a2e 



</nh 2 



At this point, we use the classical steepest descent method in order to choose 
the contour 7 such that (A.29) holds true. One has (1 — u 2 ^ n e lua£ / h = 
e n(log(i-u )+iru) Tjj-tffa r — a£ / n h e]0,o\. Thus, r > is a small parameter. 
The phase <l?(u) = log(l — u 2 ) + iru has a single nondegenerate critical point 
u c near 0, which satisfies, u c = ir/2 + C(r 3 ), and the critical value is equal to 
<J>(it c ) = -r 2 /4 + C(r 4 ). Moreover, one has $"{u c ) = -2 + 0(r 2 ). It is then 
easy to verify that one can select the contour 7 in Im(it) > r/4 connecting tt_ 
to u+, and such that on 7, one has both Re($(n)) < Re(3>(u c )) — Co|« — tt c | 2 
and |u| > Co(r + |u — u c |) for some Co > 0. We thus get 



(A.30) 



j < Ce n$(« c ) 



?• + ,s 



Then we get (A.29) from (A.30); one has n<f>(u c ) < —a 2 e 2 /8nh 2 , and since 
ry/n = ae/hy/n > clCq 1 ^ 2 , one has e _ns2 ^q^j < C'/ry/n < C%. The proof 
of Lemma 10 is complete. □ 

Next, to deduce from the L 2 estimate (A. 19) the desired L°° estimate 
(A. 15), we use the following lemma. 

Lemma 11. For given ciq,Ao, Co, there exist ai,A%,Ci,p> 0,q> such 
that for e G ]0, l],n > 1 and < h < e, the following holds true: for any func- 
tion f on M which satisfies ||/||l°° < 1? II < CqE~ 2 and 



L2(B( X0 ,e/2))<C e- A °e- a ° £ ' nh , one has 



(A.31) 
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PROOF. Let r > and xo G Co°([0, 2r [) equal to 1 on [0,r ]. Set f L = 
Xo(—h 2 A. g )f and fu = f — fL- From Lemma 8, there exists and N\,h 
such that -/h = Ei )h (l - T h ) f + M,/J an d thus from (4.14) and /i^A^I = 
2(d + 2)(l-T h ), we get 

(A.32) H/hIIloc <Ch 2 e~ 2 . 

Let $o £ Co°([0)2ro[) be equal to 1 near the support of xo- One has Xo(l — 
$o) = [we use the notation xo = Xo(~ h 2 A g ), <3?o = ^o(—h 2 A g )]. By Lemma 
3 and with |A°| defined by (3.31), we get 

Xo|A h |/ = Xo|A£|$ / - 2(d + 2) X0 A h f 

(A.33) 

-2(rf + 2) X o(T^ 2 (l-$o))/. 
Since A h G £ c 7°° and [by (A.32)] ||/i~ 2 (l - $ )/||i°° < Ce~ 2 , we get 

(A.34) II|A£|/l|U~ = ||xo|AS|*o/IU» < Ce" 2 . 

By (1.6), one has \A° h \ = -(1 + h 2 A g B)A g with 5 G Therefore, one has 

\A° h \f L = |A°|$ /i = -(1 + h 2 A g B<S> )A g f L . 

If ro is small, the operator 1 + h 2 A g B&o is invertible on L°°, and thus we 
get from (A.34) 

(A.35) I|A s /l||l~ <Ce~ 2 . 

Let ip(x) G [0,1] be a smooth function with support in the ball B(xo,e/3) 
with ip{x) equal to 1 in the ball B(xQ,e/4), and such that ||V a Vl|i/°° — 
C a e~\ a \. Set F(x) = if](x)f L (x) = iP(x)(f(x)-f H (x)). Using (A.32), A g F = 
ipA g f L + [A g ,i>]f L and (A.35), we get 

\\F\\ L 2 < c(e- A ° e - ao£2 / nh2 + h 2 £ - 2+d / 2 ), 

(A.36) 

\\A g F\\ Lao <Ce~ 2 , H^IIloo < C7. 

We now conclude that (A. 31) holds true using (A.32), (A.36) and the clas- 
sical interpolation inequality, with 9 > 

(A.37) HFHioc <C\\(l-Ag)F\\ e LOO \\F\fe 6 . 

The proof of Lemma 11 is complete. □ 

By the last inequality in (A. 18) and (A. 19), the function / = M^(ip XOt£ ) 
satisfies < 1 and ||/||L 2 (B(a;o,e/2)) < Ce~ a£2 / nh ' 2 . Let us show that it 

satisfies also |||A/j|/||£oo < Ce~ 2 . Let us recall that the operator |A/J is 
defined in (4.7). By (4.1) and (4.5), one has \ A h \ = \A h \ +2(d + 2)h~ 2 R h and 
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ll-Rftlli 00 < Ch 3 . One gets easily from (2.17) || | Ah\<Px ,e\\L°° < Ce 2 . Thus, 
one has also HIA/J^a,,, e ||£oo <C(e~ 2 + h) <C'e~ 2 . Since jA^J commutes with 
M h , one has M%(\A h \(p XOj£ ) = \A h \M%(ip XOtE ), and this implies since M h is 
Markovian, |||A\|M£(^ , £ )|| L oo < Ce" 2 . Thus we get |||A /l M^(^ 0i£ )|| L ^ < 
C(e~ 2 + h) < Ce- 2 . Prom Lemma 11, (A.16), (A. 18) and (A.20) we thus 
get, for some a, A,p,q > 0, 

V X0 , h (d g (X2 ,x ) >e)< C{e- A e- ae2 ' nh2 + hPe~% 

(A.38) 

V X0)h (d g (X2 ,x )>e) < Ch~ A e-^l nh \ 

Let a be such that < a < a/ A. It remains to observe that (A.38) implies 
(A. 15), using the second line in case h > e~ a£ ' nh and the first one if h < 

2 j 1*2 

e -ae /rin -j^g p roo f Q f Proposition 1 is complete. □ 

With the result of Proposition 1, the proof of Theorem 4 follows now 
the classical proof of weak convergence of a sequence of random walks in 
the Euclidean space M. d to the Brownian motion on M d , for which we refer 
to ([12], Chapter 2.4). Let T > be given. One has, for < 5 < c e 2 and 
he]0,h ], 

Px ,h& < 1 ^ h ~ 2T > ( l ~j) h2 < 8,d g (Xi ,X l X0 ) > 4e) 
(A.39) < % sup V y(hh (3j < I < h- 2 5,d g {Xl ,X l yo ) > As) 

y eM 

< % sup V yo , h (3j < h- 2 5,dg(Xi ,y ) > 2e) 

y €M 

2C 

<-=- sup V Zo>h (dg(X2 ,z ) >e) 

z &M,nh 2 <8 

In fact, for the first inequality in (A.39), we just use the fact that the interval 
[0, T] is a union of ~ C/5 intervals of length 5/2. The second inequality is 
obvious since the event {3j < I < hT 2 5,d g (X^ , X yo ) > Ae] is a subset of 
{3j < h~ 2 5, d g (Xy , 7/o ) > 2e}. For the third, we use the fact that the event 
A = {3j < h~ 2 5,d g (X yo ,y ) > 2e} is contained in B\J j<k (CjnDj) with B = 

{d g (X^ ,y ) > e} (k is the greatest integer < 5h~ 2 ), Cj = {d g (X yo , Xy Q ) > e}, 
Dj = {d g {Xy ,yo) > 2e and d g (X yo ,yo) < 2e for I < j}, and the fact that Cj 
and Dj are independent. 
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Using the definition (A. 4), (A. 5) of the map j Xo ,h, we get easily from 
(A. 39) the convergence for T > and e > 0, 

(A. 40) lim (lim sup -P Xn h ( max d a (u(s), u(t)) > e ) ) = 0. 

Therefore, the family of probability P xo ,h is tight, hence is compact by the 
Prohorov theorem. It remains to verify that any weak limit P XQ of a sequence 
Px ,h k , h-k — > 0, is equal to the Wiener measure W Xo . By Theorem 4.15 of [12] 
we have to show that for any m, any < t± < ■ ■ ■ < t m , and any continuous 
function f(x±, . . . ,x m ), one has 



lim / f(u(h),...,u(t m ))dP Xohk 

k—>oo J 



(A.41) = J f(xi,...,x m )p tm -t m _ 1 {x m ,x m - 1 )---pt 2 -t 1 (x2,x 1 ) 

x p tl (xi,xo)dgX 1 dgX 2 - ■ ■ d g x m . 

As in [12], we may assume m = 2. For a given t > 0, let n(t,h) G N be 
the greatest integer such that h 2 n(t,h) < (d + 2)t. By (A. 4), (A. 5), one has 

dist(w(t), Kif) < h and therefore P XOth (dist(u(t), } ) > e) = for h < 
e. Thus we are reduced to prove 

f(^ h \K i iT ) )dv X0 , h 

(A.42) 

f(xi,x 2 )pt 2 -t 1 (X2,xi)p tl (xi,x ) dgXi d g x 2 . 



From (A. 3), one has 



(A.43) 

' f(x 1 ,x 2 )Mf 2 > h) - n{tl > h \x 1 ,dx 2 )Mf 1 > h \x ,dx 1 ). 



By (A.42), (A.43), we have to show that for any continuous function f(x\,x 2 ) 
on the product space M x M, one has 



lim / f(x 1 ,x 2 )Mf^- n ^ h \x u dx 2 )Mf^ h \x ,dx 1 ) 



MxM 



(A.44) 



'MxM 

or, equivalently, 



f(xi,x 2 )pt 2 -t 1 (x 2 ,xi)p tl (xi,x ) dgXi d g x 2 , 
lhnMf^\Mf^ ) - n ^ h \f(x 1 ,-))(x 1 ))(x ) 



(A.45) 

= e t ^/\e^ A ^(f(x 1 ,-))(x 1 ))(x ). 
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Since \\M™ {t,h > \\ L oo < 1 and ||e tA 9/ 2 || L oo < 1, the following "central limit" 
theorem will conclude the proof of Theorem 4. 

Proposition 2. For all f e C°(M), and all t > 0, one has 
(A.46) lini , \\e tA ^(f) -Mf > h \f)\\ L „ = 0. 

n— >U 

Proof. Since one has \\M™ {t ' h) \\ La o < 1 and ||e* A 9/ 2 || L oo < 1, it is suffi- 
cient to prove that (A.46) holds true for / G V, with V a dense subset of the 
space C°(M), and therefore we may assume that / = ej is an eigenvector of 
A g . We set n = n(t, h), and we use the notation of Section 4.2. From (4.36) 
and n(t, h)^> 1/h, we get for some a > 0, 

(A.47) \\M n h f h \e 3 )\\ L -<Ce- at ' h2 . 
One has 

g h = (M^' h) + U )e J 

(A.48) 

= Yl (l-h 2 T k (h)) n{t ' h) e£(x) e^y)e j {y)d g y. 

T k (h)<(l-8)h-i Jm 

Let Aj = {k; \fk{h) — 2^5+2) I — £ } with e small. Then from (4.8) and Theo- 
rem 2, one has %Aj = vtij = dim Ker(A g + Aj), and for any k ^ Aj, \ J M e£(y) x 
ej(y) d g y\ < Cf.h. Using (4.9), one has \rk(h) — 2 (d+2) I — ^kh f° r an y given k. 
Take N large and split the sum in (A.48) in the two pieces Tk{h) < -/V and 
T k (h) > N. Using the L°° estimate (4.12) and the Weyl estimate (4.11) to 
bound the contribution of the sum on Tk{h) > N, we get that there exists 
C, a > and for all N, a constant C(N) such that 

(A.49) \\g h - e-^^U^e^Uoo < hC{N) + Ce~ atN , 

where ILj/j is the orthogonal projector on the vector space spanned by 
the for k G Aj. Let LT| be the orthogonal projector on Ker(A g + Aj). 
From (4.8) and Theorem 2, one has \\Hjh ~ ^j\\L 2 — Cjh. From (4.24), 
one has He^ < C(l + Tk(h)) a for some p* > d, a > 0. This implies 

\\^j,h — TEj H^^tf/Lp, < Cj, and by interpolation \\Hj t h — n^lli 2 -*-!, 00 < 
for some > 0. Then (A.49) implies 

(A.50) \\g h - e- tx ^ 2 ej\\ L co < Cjh} 1 + hC(N) + Ce~ atN . 

Clearly, (A.47) and (A.50) imply (A.46). The proof of Proposition 2 is 
complete. □ 
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